Astronomy & Astrophysics manuscript no. struve 


©ESO 2010 


September 6, 2010 





Cold gas in massive early-type galaxies: The case of NGC 1167 

C. Struve* 1 - 2 , T. Oosterloo 1 ' 2 , R. Sancisi 23 , R. Morganti 1 - 2 and B.H.C. Emonts 4 



o 

<N 

Or 

CO 



O 

U 

6 

U 
•*-> 

c3 



> 
00 

in 

O 

On 
O 
O 



1 Netherlands Institute for Radio Astronomy, Postbus 2, 7990 AA, Dwingeloo, The Netherlands 

2 Kapteyn Institute, University of Groningen, Landleven 12, 9747 AD, Groningen, The Netherlands 

3 INAF - Observatorio Astronomico di Bologna, via Ranzani 1, 40127 Bologna, Italy 

4 Australia Telescope National Facility, CSIRO Astronomy and Space Science, PO Box 76, Epping, NSW, 1710, Australia 



Received 30 June 2010; accepted 31 August 2010 



ABSTRACT 



We present a study of the morphology and kinematics of the neutral hydrogen in the gas-rich (M ffl = 1.5 X 10 10 M Q ), massive early- 
type galaxy NGC 1 167, which was observed with the Westerbork Synthesis Radio Telescope (WSRT). The H i is located in a 160 kpc 
disk (« 3 x D25) and has low surface density (< 2 M Q pc -2 ). The disk shows regular rotation for r < 65 kpc but several signs of recent 
and ongoing interaction and merging with fairly massive companions are observed. No population of cold gas clouds is observed - in 
contrast to what is found in some spiral galaxies. This suggests that currently the main mechanism bringing in cold gas to the disk is 
the accretion of fairly massive satellite galaxies, rather than the accretion of a large number of small gas clumps. NGC 1167 is located 
in a (gas-) rich environment: we detect eight companions with a total Hi mass of ~6 x 10 9 M© within a projected distance of 350 kpc. 
Deep optical images show a disrupted satellite at the northern edge of the H 1 disk. The observed rotation curve shows a prominent 
bump of about 50 km s _1 (in the plane of the disk) at r « 1.3 x R25. This feature in the rotation curve occurs at the radius where the 
H 1 surface density drops significantly and may be due to large-scale streaming motions in the disk. We suspect that both the streaming 
motions and the H 1 density distribution are the result of the interaction/accretion with the disrupted satellite. Like in other galaxies 
with wiggles and bumps in the rotation curve, H 1 scaling describes the observed rotation curve best. We suggest that interactions 
create streaming motions and features in the H 1 density distribution and that this is the reason for the success of H 1 scaling in fitting 
such rotation curves. 
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1. Introduction 

Observations show that early-type galaxies (ETGs) con- 
tinue to form their structures over cosmic times down to 
the presen t day. Evidence c omes from studie s of i onised 
gas (e.g. ISadler & Gerhardl Il985l: ISarzi et all I2006I) fine 
structure in optical morph ology (e.g. iMalin & Carted. 11983b 



ISchweizer & SeitzeR 1 19921) a n d from stellar-popula ti on in - 
vestigations (e.g. iTrageretall l2000l: iTadhunter et all [2005). 
Numerous studi es of the properti e s of Hi in and around 
ETGs (see e.g.. iKnapp et all 1 19851: iRoberts & HavnesL 1 19941: 
\ iGrossi et al ., 2009) are strongly supporting the picture of an on- 
going assembly of ETGs. About 50% of the field ETGs have 
detectable amounts of H 1, sometimes forming large, regular ro- 
tatin g disk structures (e.g. lMorganti et all 120061: [Oosterloo et all 
I2007bl) . 

In order to be able to evolve, ETGs need to keep on accreting 
cold gas. For instance, cold gas is required to explain the stel- 
lar kinematical properties obser ved in ETGs (lEmsellem et all 
l2007h . lSerra & Oosterloo! d2010b find evidence that up to 5% of 
the stellar mass in ETGs is containe d in young stars that fo rmed 
during the past ~1 Gyr. Moreover. I S err a & Oosterlool (1201 Ol) ar- 
gue that a large fraction of ETGs have accreted about 10% of 
their baryonic mass as cold gas over the past few Gyr. 

In addition, theoretical work suggests that cold gas is an es- 
sential ingredient for the evolution of ETGs. While dissipation- 
less (gas-free or "dry") mergers can explain boxy, slowly rotat- 



ing massive galaxies (e.g. iFaber et all 1 19971: iKormendy et all 
2009), a comparison of N-body simulations with observations 
suggests that present-day ETGs can not h ave assembled mor e 
than 50% o f their mass via dry merg i ng (iNipoti et all l2009h . 
iJesseit et al.1 (120071) and iHopkins et all (120091) have shown that 
simulations of dry mergers are not able to reproduce some ob- 
served quantities such as counter-rotating cores, velocity dis- 
persion features and structural properties of ETGs, whereas the 
inclusion of a gas component in the simulations (10% of the 
baryonic mass) helps to reproduce the main observed features. 
Moreover, dissipative mergers and satellite accretion events are 
neede d to explain, e.g., the dynamical structure of disk-like 
ETGs dNaabet all [2006b . 

The accretion/merging of fairly massive companion galax- 
ies is one possible process to provide the required gas 
masses to ETGs, but other mechanisms can be important as 
well. In particular the slow but long-lasting infall of cold 
gas from the IGM is predicted as a possible mechanism 
which c an provide large g as masses over cosmic times (e.g. 
iBinnevt 119771: iKeres et all 120051) . In the framework of hi- 
erarchical galaxy formation, new cold gas is added to the 
halo or deposited in the outer parts of galaxies where it 
forms a reservoir to fuel the inner parts of galaxies. However, 
the exact mechanism, how the gas cools (or remains cool), 
condenses and penetrates to the inner disk with out getting 
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shock heated is currently still under debate (e.g. iPeek etal 
20081: iDekel & Birnboiml 120081: [Kawata & Mulchaevi 12008 



Broo ks et all 120091 : IKer es & Hernquist, 2009). 
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Whether the accretion/merging of a few fairly massive satel- 
lites is the main disk building process or whether the continuous 
infall of many small gas clumps plays a dominant role for the 
continuing assembly of ETGs is not clear. Given the importance 
of cold gas for the evolution of ETGs, motivated by observa- 
tions and theoretical work, it is essential to investigate the origin 
of cold gas in ETGs. 

The question of how galaxies accrete gas is not only impor- 
tant for ETGs, but also for spiral galaxies and hence for galaxy 
evolution as a whole. Signs of gas accretion are found in a large 
number of galaxies (for instance: gas-rich minor mergers, lob- 
sided disks, gas fi laments, extra-planar Hi; for a review, see 
ISancisi et al.Ll2008l) . showing the continuing assembly of galax- 
ies. For spiral galaxies the continuous supply of cold gas is im- 
portant because gas consumption times are usually shor t in th e 
inner regions (typically < 1-2 Gyr; see, e.g. jBigiel et alll2008l) . 
Hence, these galaxies can only maintain their observed star for- 
mation rates if they manage to keep on acquiring gas, either 
through merging or via the accretion of small gas clumps from 
the IGM. An interesting aspect in that respect is whether ETGs 
accrete gas in the same way, or if other mechanisms might sup- 
press the accretion of cold gas onto the disk, which may then 
explain the absence of current star formation. 

Neutral hydrogen studies of gas-rich, massive ETGs provide 
a powerful tool to investigate where the cold gas in these sys- 
tems comes from. Because the cold gas column densities must be 
low (since large-scale star formation does not occur), the H i in 
H i-rich ETGs i s often distributed in large str uctures (see, e.g., 
iMorganti et"all 120061: lOosterloo etail l20Q7bh . At large radii, 
dynamical timescales are large, allowing to investigate the evo- 
lution of the galaxy over the past few Gyr. Since ETGs con- 
tinue to form their structures down to the present day, the role 
of the H i in this formation process can be investigated over Gyr 
timescales. 

A second reason to observe H i-rich ETGs is that it allows 
to study in detail the accretion of small gas clumps. In recent 
years, neutral hydroge n gas has been dete cted in the halo of spi- 
rals (for a review, see ISancisi etaL l l2008h and the existence of 
H i in the halo could be partly explained by the accretion of small 
cold gas clumps. However, an alternative mechanism for bring- 
ing H i in the halo is star formation (supernova explosions and 
stellar winds blow gas into the halo) and it is not easy in spiral 
galaxies to determine which of the two mechanism is the most 
important one. 

To study the accretion of small gas clouds, it is thus use- 
ful to observe galaxies with no or little star formation, such as 
early-type and low surface brightness galaxies that have no/few 
supernovae e xplosions. Only very few such galaxies have been 
studied so far. iMatthews & Woodl(l2003l) detected H i in the halo 
of the LSB galaxy UGC 7321 but concluded that even in this 
galaxy the energy provided by supernovae explosions might be 
sufficient to explain the gas in the halo. 

In this paper we investigate how the cold gas in the early-type 
(SAO) galaxy NGC 1 167 (UGC 2487) is accreted in order to un- 
derstand where the H i comes from. We present deep H i observa- 
tions performed with the Westerbork Synthesis Radio Telescope 
(WSRT). NGC 1167 is a giant disk galaxy (Af B = -21.7 mag, 
D25 = 56 kpc), seen at low inclination, known to host a relatively 
powerful (10 24 W Hz" 1 at 1.4 GHz) radio source (B2 0258+35). 
What makes this galaxy interesting for the study of the origin 
of the cold gas is that previous Hi observations had shown a 
regular rotating 160 kpc d i sk with a total Hi m ass > 10 10 M Q 
dNoordermeer et all 120051: lEmonts et all 12007b . The Hi sur- 
face density is low (< 2 M© pc -2 ) and the Hi disk extends 



Table 1: Optical and radio parameters for NGC 1 167. 



Parameter 


Value 


Reference 


Morphological type 


SAO 


1 


Center a (J2000.0) 


03 01 42.4 


1 


Center 6 (J2000.0) 


35 12 20.7 


1 


Distance (Mpc) 


70 


1 


L B (L q ) 


7.4 x 10 iU 


2 


Systemic velocity (km s _1 ) 


4951 +3 


3 


D 25 (kpc) 


49 


1 


An (kpc) 


160 


3 


Total H 1 mass (M©) 


1.5 x 10 10 


3 


H 1 inclination (deg) 


38 


3 


Mean P.A. (no warp) (deg) 


75 


3 


Total mass (M Q ) 


> 1.6 x 10 12 


3 



Note. - Units of right ascension are hours, minutes, and seconds 
and units of declination are degrees, arcminutes, and arcseconds. 
References. - (1) NED; (2) based on LEDA; (3) this work. 



Table 2: Summary of observations. 



Total observing time 


138 h 


Total on- source integration time 


~ 107 h 


Observing period 


Oct. 2006 to 




Nov. 2008 


Bandwidth 


20 MHz 


Number of spectral channels 


1024 


Channel separation 


8.5 km s _1 


Velocity resolution (after Hanning) 


17.0 km s- 1 


Bandpass/flux calibrator 


3C48 and 3C147 


Weighting scheme 


robust 0.4 


Resolution 


28'.'9 x 16V8 


Beam position angle 


\°.5 


rms noise (mJy beam -1 ) in the cubes 


0.125 


rms noise per channel (10 18 atoms cm -2 ) 


4.8 


rms noise per channel (M© pc ~ 2 ) 


0.039 


Detection limit (3<x over 3 channels) (10 6 M o ) 


5.6 



beyond the optical disk out to ten 7?-band disk scale lengths. 
Deep optical observations show a faint , tight ly wound spiral 
structure at r < 30 kpc (lEmonts et all l2Q10h . No stellar bar 
structure is visible in the data. Ha emission-line imaging re- 
vealed no emission (log L < 37.6), except for close to the 
AGN (r < 5 kpc), and he nce no sign of ongoing star formation 
(Po gg e~& EskridgeU 19931) - in agreem ent with Ha long-slit ob- 
servations obtainedbyjEmon^ observations at 2 
-10 keV (lAkvlas & Geor gantopo ulosl 120091) show a X-ray halo 
which extends beyond D25 (J. Rasmussen, private communica- 
tion). In addition to the investigation of the origin of the cold gas, 
the very extended H 1 disk allows to derive a rotation curve out 
to large radii and to investigate the mass distribution out to large 
distances from the centre. Some optical and radio parameters are 
given in Table [T] 

The organisation of this paper is as follows: in Sect. 2 we de- 
scribe the observations, while Sect. 3 presents the observational 
results and in Sect. 4 we search for a possible H 1 halo compo- 
nent. In Sect. 5 we derive a rotation curve to investigate the disk 
dynamics. A discussion is given in Sect. 6 and a summary in 
Sect. 7. The redshift of NGC 1167 (z = 0.0165) translates to a 
distance of D = 70 Mpc assuming Hq = 71 km s" 1 and conse- 
quently 3" correspond to about 1 kpc. 
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Fig. 1 : DSS2 image showing the H i group environment of NGC 1 167. Numbers correspond to numbers in Table[3] H i contour levels 
are: 1, 2, 4, 8, 16 and 32 x 10 19 atoms cm -2 . 



2. Observations and data reduction 

A ten-hour WSRT service observation in 2006 was followed-up 
by eleven 12-h observing runs in 2008 (see Tab. [2]). The net on- 
source integration time is about 107 h. The strong radio contin- 
uum source in the centre (1.8 Jy) required non-standard calibra- 
tion. In order to guarantee a stable bandpass, calibration observa- 
tions of 15 minutes were performed every two hours. 3C48 and 
3C147 were used as calibrators. An almost uniformly filled uv- 
plane could be achieved for the combined data by interleaving 
the hour angle of the calibration observations between individ- 
ual observations. 

The observing band was centred on the redshift of 
NGC 1167. The total bandwidth was 20 MHz, corresponding 
to a velocity range of -3700 km s" 1 , and was covered with 1024 
channels (dual polarisation). In order to increase the signal-to- 



noise ratio, two channels were binned and subsequently Hanning 
smoothed resulting in a spectral resolution of 17 km s" 1 in the 
output data cube. The data reductio n was performed using the 
MIRIAD package dSault et al.L[l995h . 



After flagging, bandpass and phase calibration were per- 
formed. In addition, on every individual data set self-calibration 
of the continuum emission was done and the results of this also 
applied to the line data. The continuum model was subtracted 
from the line data and the residual continuum emission was re- 
moved by making a fit to the line-free channels of each visibility 
record. The data cube was cleaned using the Clark algorithm. For 
the extended emission a cube with robust weighting of 0.4 gives 
the best compromise between spatial resolution and sensitivity, 
which is used for the analysis of the H i data. 
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Table 3: Properties of H i detections. 



No. 


Identification 


RA 


Dec 


^sys 


Av 


Av off 


D 


A>ff 








M m 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


(13) 





NGC 1167 


03 01 42.4 


35 12 20.7 


4951 


520 


- 


69.7 


- 


- 


457.65 


461.33 


15.324 


1 


J0302478+3504284 1 


03 02 47.8 


35 04 28.4 


4864 


103 


-87 


68.5 


15.5 


310.4 


1.62 


2.54 


0.081 


2 


2MASX 2 


03 01 48.6 


35 05 42.4 


4903 


213 


-48 


69.1 


6.7 


134.8 


16.29 


17.70 


0.578 


3 


UGC 2465 


03 00 37.4 


35 10 08.7 


5065 


529 


114 


71.3 


13.4 


267.6 


57.23 


80.90 


2.812 


4 


J0302091+3511120 3 


03 02 09.1 


35 11 12.0 


4907 


34 


-44 


69.1 


5.6 


111.6 


0.39 


0.41 


0.013 


5 


J0301222+3514055 4 


03 01 22.2 


35 14 05.5 


5057 


42 


106 


71.2 


4.5 


89.4 


0.78 


0.81 


0.028 


6 


J0302095+3516308 3 


03 02 09.5 


35 16 30.8 


4903 


77 


-48 


69.1 


6.9 


138.0 


13.20 


14.39 


0.470 


7 


J0302307+3516581 4 


03 02 30.7 


35 16 58.1 


5245 


128 


194 


73.9 


10.9 


217.9 


6.07 


7.56 


0.282 


8 


J0301491+3529012 1 


03 01 49.1 


35 29 01.2 


4954 


94 


3 


69.8 


16.7 


334.0 


25.59 


43.10 


1.436 



Note. - (1) Number in Fig.Q] (2) Name (NED)/Identification. (3) Centre: Right ascension. (4) Centre: Declination. (5) Systemic ve- 
locity [km s" 1 ]. (6) Line width above 3<x. (7) Velocity off-set from NGC 1 167. (8) Distance from Earth assuming Hubble flow (Hq = 
71) [Mpc]. (9) Off-set from NGC 1 167 [arcmin]. (10) Off-set from NGC 1 167 [kpc]. (11) Integrated flux [Jy km s" 1 ]. (12) Primary 
beam corrected Flux [Jy km s -1 ]. (13) H i mass [10 9 M o ]. 1 counterpart on DSS image; 2 full name is: 2MASX J03014864+3505423; 
3 star in front; 4 no counterpart on DSS image. 
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Fig. 2: Radial H i surface density profile. 
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Fig. 3: 5-band image from lEmonts et al.l (1201 Oh showing the 
faint spiral structure of the stellar disk. 

3. Results 

3.1. NGC 1167 

The high sensitivity of our observations has allowed to de- 
tect new faint Hi structures compared to earlier observations. 



Figure \T\ shows the total Hi image of NGC 1167 and its neigh- 
bours. This image was obtained by smoothing the data cube to 
35" resolution and using this smoothed cube to create a mask 
(blanking all signal below +3<x rms in the smoothed cube) which 
was applied to the original cube. The Hi extends over 160 kpc 
in diameter (« 3 x D25), with a total Hi mass detected in emis- 
sion of 1.5 x 10 10 M Q (which is roughly thre e times the H 1 mass 
of the Milky Way; Henderso n et all Il982l) . Channel maps are 
shown in Fig. I A. II and the radial Hi density profile in Figure[2 
The highest H 1 surface density is found at the radius of the faint, 
tightly wound spiral structure found in optical images (Fig. [3j 
see also Emonts et al. 2010). Absorption is detected against the 
strong central continuum, which is unresolved at our resolution. 
Higher resolution observations show that the absorption is de- 
tected against a kpc scale structure. This will be discussed in a 
future paper (Struve et al. in prep.). The apparent central drop in 
surface density (Fig. [2]) is affected by the absorption. Despite the 
large H 1 mass, the surface density is too low for star formation 
to occur globally. 

To the southwest of NGC 1 167, the H 1 extends to larger radii 
as compared to the rest of the galaxy (see Fig. |4l left panel). 
The gas in this region does not follow exactly the rotation of the 
regularly rotating ~65 kpc disk (see also Fig. lA.lb . 

The derivation of a velocity field is not straightforward. 
Some of the velocity profiles are broad (see Fig. [5] and [7] for a 
few examples) and close to the centre (r < 100") beamsmearing 
becomes severe due to the large beam size (9.6 x 5.6 kpc) and the 
steep velocity gradients. To minimise beamsmearing effects, we 
use the velocity of the peak of the profile instead of the intensity 
weighted mean. The resulting velocity field (Fig.|4l right panel) 
shows that the gas kinematics for r < 65 kpc is dominated by 
regular large-scale rotation. 

The radial velocities along the major axis and the main rota- 
tional properties are shown in Fig. The rotation curve (derived 
in Sect. 5 and over-plotted on the major- axis slice in Fig. [5]) is 
rising within the central beam and is essentially flat out to the 
largest radii measured. However, on both the approaching and 
receding side, at radii between ~26 and -45 kpc, the rotation 
curve makes a 25 km s" 1 jump up and goes down again by about 
37 km s" 1 (uncorrected for inclination). The radius of the peak of 
the jump (r « 36 kpc) and the radial extent of this jump are sim- 
ilar, but not completely identical for the approaching and reced- 
ing side. This "bump" in velocity occurs just outside the bright 
optical disk at r « 1.3 x R25, exactly at the radius where the 
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Fig. 4: Left panel: Total Hi intensity map of NGC 1167. Contour levels are: 1, 2, 4, 8, 16 and 32 x 10 19 atoms cm" 2 . Right panel: 
Velocity field based on the peak of the profile (see text). The white contour indicates the systemic velocity, 4951 km s" 1 , the black 
and grey contours give the velocity in increments of 50 km s" 1 . The approaching side is to the East. The beam is indicated in the 
lower left corners. 
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Fig. 5: Position- velocity slice along the major axis. The white squares give the rotation curve as derived from modelling (see 
Sect. 5). The bump in the rotation curve is clearly visible. Contour levels: -9, -6, -3, -2 (thin black), 2, 3, 6 and 9cr (thick black). 
The horizontal dashed line indicates the systemic velocity. 



Hi surface brightness rapidly drops (Fig. [2]). This is further dis- 
cussed in Sect. 6.3. 

3.2. The (H \) environment of NGC 1 167 

We detect eight sources (Fig. [T]) in the vicinity of NGC 1167 
with Hi masses between ~ 1 x 10 7 M© and 3 x 10 9 M© within a 
projected radius of 350 kpc. The total amount of H i in the envi- 
ronment is 5.7x 10 9 Mq. Two of the detections are identified with 
catalogued galaxies. Table [3] summarizes the basic properties of 
all Hi detections. 



The most massive galaxy besides NGC 1167, UGC 2465, is 
located 268 kpc (projected) west and shows a distorted H i mor- 
phology and kine matics. In addition , deep optical observations 
(for details, see lEmonts et all 120101) show a disrupted satellite 
at the northern edge of the H i disk with a possible tail towards 
the region of the southwestern H i extension (Fig. [6j see also 
Sect. 6.1). Besides these two cases, none of the other Hi detec- 
tions show hints of an interaction with NGC 1167 in terms of 
bridges or tails. 

Five detected galaxies (No. 1, 2, 6, 7 and 8 in Tab. 0) show 
indications of regular rotation, despite the low spatial resolution. 
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Right Ascension (J2000) 

Fig. 6: #-band image from iEmonts et al.l (1201 Ol) overlaid 
with Hi density contours. Contour levels: 1, 8 and 16 x 
10 19 atoms cm -2 . The arrow indicates the location of the dis- 
rupted satellite. 

Objects No. 2, 6 and 7 are fairly gas-rich with Hi masses above 
10 8 M Q . DSS2 images do not show an optical counterpart for 
No. 6 and 7 (see Fig.Q]). However, the DSS2 images are shallow 
and optical counterparts with luminosities several times 10 7 L© 
could still be present. Close to NGC 1167, at distances of 90 to 
110 kpc from the centre, we detect two unresolved Hi clouds 
(No. 4 and 5) with masses just above our detection limit which 
do not have an optical counterpart in DSS2 either. No additional 
low-mass H i clouds were found i n a sy stematic search using the 
source finder Duchamp ( Whiting J2008I) . 

4. The search for an H i halo 

One of the main reasons to observe NGC 1167 is to investigate 
the origin of the large amount of H i in this early-type galaxy. 
One possibility is the long-term accretion of a large number 
of small cold gas clumps. Signs for the accretion of such gas 
clouds have been found in the ha los of a few spiral galaxies (for 
a review, see lSancisi et all 120081) . Over a Hubble time, galaxies 
could accrete a significant fraction (> 10 9 M ) of their cold gas 
in this way. In order to investigate the possibility that NGC 1 167 
also accretes such gas condensations, we search for an extended 
halo of Hi clouds. In Sect. 4.1 we explain how to identify gas 
that is located in the halo and in Sect. 4.2 we show that no 
H i clouds are detected in the halo of NGC 1 167 and place upper 
limits to the accretion of small gas clumps. 

4.1. Identifying halo gas 

In low-inclination galaxies, such as NGC 1167, the separation 
of gas located in the disk and the halo is not straightforward. 
This is because gas located in the halo lies along the same 
line-of-sight as disk gas, unlike in edge-on systems. However, 
the kinematics of the halo gas differs from the disk kinemat- 
ics such that the velocities of the gas can be used to separate 
the halo component from disk gas. For instance, it has been 
found that gaseous hal os observed in spirals are lagging with re- 
spect to disk rotation dFraternali et all 120021: Ikeald et all 120071: 



lOosterloo et all l2007al) . This lagging gas is most clearly seen 
in a /?v-diagram along the major axis in the form of wings on 
the side of the lower rotation velocities, toward systemic. These 
systematic asym metries in the H i velo city profiles have been ref- 
ered to as beards dSancisi et alll200l1) . Low-inc lination galaxies 
that s how prominent b eards are e.g. NGC 240 3 dFraternali et all 
120021) and NGC 6946 (iBoomsma et alll2008l) . 

Position- velocity slices along the major axis of NGC 1167 
(Fig. [5} or near the major axis (Fig. [3 right panel) show H i line 
profiles that are strongly asymmetric with respect to the peak 
of the line, in particular in the inner part of the disk (similar to 
NGC 2403 and NGC 6946). Note that the white squares, which 
mark the rotation curve (which we derive in Sect. 5), follow 
closely the peaks of the profiles. In some locations the profiles 
extend (as a beard) 200 km s" 1 toward the systemic velocity 
(Fig. 13 right panel). 

However, there is a major difference between NGC 1 167 and 
previously studied low-inclination galaxies: NGC 1167 is lo- 
cated relatively far away such that the linear resolution of our 
observations is coarse (beam is 9.6 x 5.6 kpc). As a consequence 
the velocity gradients within a beam are large. Hence, broad 
velocity profiles are observed, in particular close to the centre 
where the velocity gradients are largest. Because of the large ro- 
tation amplitude, the velocity gradients would even be larger if 
NGC 1 167 would be observed at the same resolution as the pre- 
viously studied objects. The convolution of the signal with the 
observed beam spreads the emission to neighbouring locations 
such that the broad profiles observed in Fig [5] and [7] could be the 
result of beamsmearing instead of indicating a halo. To identify 
a possible H i halo component, it is therefore important to com- 
pare the data to a model that takes beamsmearing into account. 
In addition to the velocity gradients near the centre, NGC 1167 
shows evidence for large, local velocity gradients at other loca- 
tions due to streaming motions in the disk (see Sect. 5). Also 
these require that beamsmearing is taken into account in mod- 
elling the H i disk of NGC 1 167. 

The thin disk of NGC 1167 is modeled by placing at every 
pixel in the cube a Gaussian profile centred at the velocity at that 
location as derived from the observed velocity field. This cube 
is then convolved with the beam. Thus, this approach ensures 
that the local velocity structure and beamsmearing effects are 
taken into account. The intensity of the profile is taken from the 
total H i map at the position. For the Gaussian profile, a velocity 
dispersion of 9 km s" 1 is assumed. 

4.2. An upper limit to the H i halo 

Comparing the model cube with the data shows that the model 
describes the overall kinematics of the H i in NGC 1 167, includ- 
ing the broad profiles seen in many locations. This is illustrated 
in Fig. |7] where we show two representative pv slices. The fact 
that our beamsmeared model of the thin H i disk does explain the 
broad profiles indicates that these are due to beamsmearing and 
that we do not observe Hi clouds in the halo of NGC 1167, in 
contrast to what is seen in many later- type spiral galaxies. 

Figure |3 also illustrates that the only features not consistent 
with our model are at the noise level. However, many other such 
"blobs" are present in the data cube, also at locations and veloci- 
ties quite remote from NGC 1 167. Therefore, we take the masses 
and column densities of the clouds of the kind indicated in Fig. [3 
as the upper limit for any halo population in NGC 1 167. 

The upper limit of the peak surface density of a possible 
Hi halo component is hence 0.05 M© pc" 2 , which corresponds 
to a column density of 6 x 10 18 cm" 2 and any cloud in the halo 
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Fig. 7: Two east- west /?v-slices of the data cube overplotted with the l<x contour level (red contour) of our model. Contour levels of 
the data: -3, -1.5 (grey), 1.5, 3, 4.5 and 6<x (black). H i emission not covered by the model is indicated by arrows. 



of NGC 1167 has an Hi mass of less than 10 7 M Q . We note, 
for comparison, that in some spiral galaxies halos have been 
observed with average column densities well above 10 19 cm" 2 
while the regions with highest Hi surface density in the halo 
of, e.g., N GC 2403 and NGC 6946 have column densitie s above 
10 20 cm" 2 dFraternali et al.l.l2002l:lBoomsma et al.Ll2QQ8h . which 
is a factor ~20 above the upper limit for NGC 1167. We can, 
therefore, rule out that NGC 1167 has an Hi halo with similar 
(column) densities as compared to those late-type, star-forming 
galaxies. The mean density contrast ratio between H i located in 
the disk and in the halo is > 20 for NGC 1167. This is signifi- 
cantly larger than for NGC 2403 and NGC 891 (ratio is 10 and 
2.3 respectively), but is similar to the one of NGC 6946 (ratio is 
22). 

5. Rotation curve 

In this section we further investigate the nature of the veloc- 
ity bump seen in the kinematics of NGC 1167, as described 
in Sect. 3.1. In particular, we investigate whether this fea- 
ture reflects the overall mass distribution of this galaxy, or 
whether it could be indicative of a recent disturbance suffered 
by NGC 1 167. We derive a rotation curve and show that none of 
our applied mass models explains the feature in the kinematics 
satisfactorily. We conclude that large-scale streaming motions 
must be present in the disk. 

The rotation curve is derived by fitting a set of tilted rings 
to the observed velocity field (Fig. ffl right panel) (following 
Beg emanll 19871) . The resulting rotation curve is plotted in Fig. [8] 
and given in Tab. IB. 11 The rotation curve has two sources of un- 
certainty. Every data point has an individual error of ~13 km s" 1 
(estimated from the residuals in the velocity field, Fig.0. In ad- 
dition, a systematic uncertainty exists wh i ch dep ends on the in- 
clination of the disk. Folio win g IB e gemanl £1987), the systematic 
uncertainty in inclination for NGC 1 167 is 2°, which is constant 
with radius (i = 3 8° ±2°). This translates into a possible system- 
atic offset of 15 km s" 1 . The position angle is constant within 
the uncertainties (PA = 75°±3°). The remarkable bump seen in 
Fig. [5] is also visible in the rotation curve as a distinct feature 
at radii between 80" and 130". In principle, the feature in the 
kinematics could also be explained as the result of a change in 
inclination instead of rotation velocity. However, a bending up 
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Fig. 8: Rotation curve. The error bar indicates the uncertainties. 



and down in inclination of ~6° would be required. Such a bend- 
ing is not observed in other galaxies and therefore we think that 
it is unlikely that this feature reflects such a change in inclina- 
tion. 



5.1. Large-scale streaming motions 

The presence of non-circular motions is also evident from the 
residual velocity field. In Figure [9] (top panel) the velocity field 
of an axisymmetric model, which is based on our derived rota- 
tion curve, has been subtracted from the observed velocity field. 
A large-scale residual structure is seen at the radius of the fea- 
ture in the rotation curve. We note that the residuals are not com- 
pletely axisymmetric, in agreement with the slightly asymmetric 
appearance of the feature in the rotation curve (Sect. 3.1). This 
residual structure is even more prominent assuming a purely flat 
model rotation curve (i.e. without a bump), as shown in the bot- 
tom panel of Fig. [2 Additional velocity deviations from the ax- 
isymmetric model at larger and smaller radii are also present. 

The systematic residuals imply that streaming motions must 
be present in the disk. Simple radial motions are excluded be- 
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Table 4: Parameters for mass models with isothermal halos, Hi scaling and MOND. 



model 


n 


T d 


1 


r c 


Po 


A 


x 1 




Mq/Lq 


M /L 




kpc 


M pc" 3 






(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


minimal^ 2 


4.1 


3.6 


1.4 


10.6 


14.5 




12.2 


max. allowed disk 


4.3 


4.7 


1.4 


25.6 


3.0 




16.3 


min. allowed disk 


3.8 


1.0 


1.4 


2.7 


249.0 




18.7 


constant M/L 


4.2 


4.2 


1.4 


21.0 


4.4 




16.6 


H i scaling 


4.5 


4.9 


46.4 








11.2 


MOND 


3.8 


4.9 


1.4 






3565 


13.8 



Note. - (1) model; (2) bulge 7?-band mass-to-light ratio; (3) disk 7?-band mass-to-light ratio; (4) Hi scaling factor; (5) halo core 
radius; (6) halo central density; (7) Mond factor; (8)^ 2 value from fit. 
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Fig. 9: Top panel: Residual velocity field based on the derived 
rotation curve (Fig. [8]). Bottom panel: Residual velocity field as- 
suming a flat rotation curve (see text). The beam is indicated in 
the lower left corners. 



5.2. Mass models 

We model the mass distribution in order to investigate whether 
the feature in the rotation curve can be explained by one of the 
standard models. We assume that the gravitational potential of 
a disk galaxy can be written as the sum of the individual com- 
ponents: stellar bulge, stellar disk, gaseous disk (multiplied with 
the H i scaling factor rj) and a dark matter halo. To model the 
stellar mass distribut ion we use published 7?-band photometry 
(lNoordermeer[ [2006) and adopt Noordermeer's disk/bulge sep- 
aration. The uncertainty of the rotation curve affects the M/L 
estimate by < 0.3. 

None of the tested models reproduces the feature in the 
rotation curve satisfactoril y (see Fig. [TUl and T ab. |4]). A trun- 
cated optical disk (see, e.g.. Ivan der Kruit,2007|) is not observed 
and, therefore, neither explains the feature in the rotation curve. 
Moreover, even if a truncation in the optical disk would have 
been observed at the right radius, it would not have been able, 
given the faintness of the disk, to explain the amplitude of the 
bump. 

The best model (smallest overall x 2 ) is that of a scaled-up 
version of the H i mass distribution with no dark matter halo. The 
H i scaling factor o f 33 (correcting for h elium) is higher than the 
average value that lHoekstra et al.l d2001) found for a sample of 
24 spiral galaxies (77 = 6.5), but is in agreement with the sys- 
tematic trend that more luminous systems require higher scaling 
factors than low luminosity systems. High scaling factors have 
also been f ound for other lumin ous galaxies, e.g., 77 = 37 for 
UGC 6787 (iNoordermeerl [2006h . However, also in the Hi scal- 
ing model the feature in the rotation curve is not perfectly repro- 
duced. The predicted bump is approximately 5 kpc closer to the 
centre than the observed bump. 

The fact that none of the tested models fits the observed ro- 
tation curve suggests that the feature in the rotation curve does 
not reflect an equilibrium mass distribution of NGC 1167. In 
Sect. 6.3 we will argue that the bump in velocity is consistent 
with the presence of large-scale streaming motions in the disk, 
as already suggested from the residual velocity field (Fig. [9] and 
Sect. 5.1). 



cause the feature in the rotation curve is seen strongest along the 
major axis (where radial motions should be zero) and not along 
the minor axis (where they should be maximal). The measured 
peak velocity of the large-scale residuals is above 34 km s" 1 
corresponding to a deprojected amplitude (i.e. velocity in the 
plane of the disk) of 55 km s" 1 . In Sect. 6.3 we argue that 
these large-scale streaming motions may be caused by an inter- 
action/merging event. 



6. Discussion 

We start this section by discussing the formation of the H 1 disk 
and the role of different accretion mechanisms (Sect. 6.1). In 
Sect. 6.2 we give possible explanations for the absence of a pop- 
ulation of cold gas clouds and Sect. 6.3 briefly comments on a 
possible connection between the H 1 density and the large-scale 
streaming motions in the disk. 
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Fig. 10: Mass models. The fitted individual contributions to the rotation curve are stellar bulge (short dashed), stellar disk (long 
dashed), gas disk (continuous line) and dark halo (dotted). The dot-dashed line shows the resulting model rotation curve. The error 
bar indicates the uncertainty for an individual data point of the observed rotation curve. 



6.1. Formation of the H i disk 

Where does the H i in NGC 1 167 come from? The various signs 
of merging and interaction activities suggest that currently the 
Hi is brought in by fairly massive satellite galaxies. This means 
that most of the H i arrives in large quantities in a few events. 

Signs of recent and/or ongoing interactions with the environ- 
ment of NGC 1 167 are, for instance, the southwestern extension 
of the H i disk, the feature in the rotation curve, the disturbed 
H i morphology and kinematics of UGC 2465 and the disrupted 
satellite at the northern edge of the Hi disk (Sect. 3.2). All of 



these signs are pointing towards an ongoing assembly of the 
H i disk. About 30% of the detected H i in the data cube is lo- 
cated in companions. Hence, it might be possible that, over the 
course of several Gyr, NGC 1 167 has built-up a large fraction of 
its H i disk via minor mergers and interactions with the environ- 
ment. 

The southwestern extension could be the result of such a mi- 
nor merger. Figure [6| shows the disrupted satellite at the northern 
edge of the H i disk with a possible tail-like structure all the way 
to the H i extension (Sect. 3.2). Therefore, we might see also the 
disrupted H i in the aftermath of a minor merger. The H i mass 
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of the southwestern extension is at least 9 x 10 8 M Q , showing 
that the accretion of satellite galaxies can bring in large amounts 
of H i. This minor merger may have also induced the large-scale 
streaming motions in the disk that are seen as the feature in the 
rotation curve (Sect. 5.1, see also Sect. 6.3). The non-circular 
motions provide, therefore, another piece of evidence for the ac- 
cretion of, fairly massive, companion galaxies. 

Another possibility is that the H i extension could be caused 
by a recent interaction with UGC 2465. In this case, the past 
interaction time estimate is roughly 1 Gyr, given the projected 
distance between both galaxies and assuming a relative velocity 
offset of 200 km s" 1 . This time estimate equals one revolution 
of NGC 1167 (at r = 65 kpc), but is less than the 1.3 Gyr that 
one revolution of UGC 2465 takes. After an interaction it takes 
a few rotations before the gas settles in a rotating disk structure. 
Hence, it is not a surprise that the southwestern extension and 
part of the H i structure of UGC 2465 is not completely settled. 

Other mechanisms that could contribute to the H i disk build- 
ing process do not seem to play a major role. Below, we briefly 
discuss the relevance of small gas clouds, the possibility of a 
major merger and the role of stellar mass loss. 

Accretion of small gas clumps: 

We do not find a large number of small gas clumps. In fact, only 
two clouds just outside the H i disk (No. 4 and 5 in Tab. 0) are 
detected. The very few tentative H i structures (with masses of 
~10 7 M Q ), which are not explained by the model (Sect. 4), are 
features at the noise level and form an upper limit for a possi- 
ble halo cloud population. The two detected clouds have masses 
at the high mass end of the high velocity clouds (HVCs) ob- 
served around the Milky Way and other nearby galaxies (10 4 <; 
M HT <10 7 Mp: e.g.lThilkeretall l2004t IWestmeier etaTL 12005b 
Wa kker et all I2007L 120081) . Clouds wit h Hi masses ~ 1 7 M Q 
were also detected around NGC 2403 dFrate rnali et al, 2002), 
NGC 891 (lOosterloo et~all 12007 ah N GC 6946 dBoomsma et all 
2008]) and NGC 2997 (|Hess_eiall 120091) . What makes these 
HVCs interesting is that they could be, at least in some cases, gas 
condensations from the IGM. The presence of two such clouds 
around NGC 1167 might, therefore, indicate that the accretion 
of small gas clouds, at a low level, occurs. 

Nevertheless, the absence of a large number of such clouds 
shows that the "diffuse" accretion of small cold gas clumps from 
the IGM is not the dominant disk building process. The number 
of required clumps (> 10 3 ) would be unphysically large, given 
that we detect only very few such gas clouds. Possible explana- 
tions for the absence of gas clouds located in the halo are dis- 
cussed in Sect. 6.2. 

Major merger: 

While gas-rich major mergers can form giant low column den- 
sity Hi disks (Sect. 1), we have no evidence in our Hi data that 
this is actually the case for NGC 1167. Neither do the deep op- 
tical images show any indication of a major merger. Such a dra- 
matic event would certainly need to be older than a few Gyr, 
given that the Hi disk shows regular rotation for r < 65 kpc. 
After a major merger, it takes at least a few revolutions (corre- 
sponding to a few Gyr at r = 65 kpc) before the gas s ettles in 
a reg ular disk structure. A stellar population analysis of lEmontsI 
( 200Q) also suggests that the last interaction (that caused a burst 
of star formation ) occured at least one to several Gyr ago, but 
this starburst may be the result of a smaller accretion event and 
is no evidence for a major merger. 



Stellar mass loss: 

The Hi can not originate from stellar mass loss. A H stars shed 
mass over their lifetimes. According to models of ICiotti et al.l 
(1199 ll) and lBrighenti & Mathewsl (l2000h . the expected cold gas 
mass from stellar mass loss processes for massive galaxies could 
be > 10 10 M Q over a Hubble time. In the case of NGC 1167 it 
is unlikely that the H i results from stellar mass loss because the 
H i disk extends well beyond the optical disk (£>hi = 3 XR25) and 
carries therefore more angular momentum than the stars. Also in 
other galaxies there is little evidence that cold gas comes from 
mass loss processes. In the cases where gas could be associate d 
with stellar mass loss, it is ionised (see, e.g JSarzi et all 2006). 

As a final remark, we would like to note that all eight detected 
companions of NGC 1 167 have systemic velocities close to that 
of NGC 1167 (none is deviating by more than 200 km s" 1 , see 
Tab. [3]). Furthermore, all but No. 7 (in Tab. 0) have systemic ve- 
locities which follow the sense of rotation of NGC 1 167. We can- 
not rule out that this i s the result of projection effects. However, 
lErickson et al.l (1 19991) found the trend that companions seem to 
be close in velocity to the host galaxy in a number of (late-type) 
spiral galaxies. They argue that their sample is not affected by 
projection effects and they find evidence for varying sizes of 
dark matter halos. However, for galaxy sizes above 200 kpc their 
constraints become weak. Larger statistical samples are needed 
to shed more light on this issue. 

6.2. The absence of an H 1 halo 

Different theoretical explanations have been proposed that could 
explain the absence of an extended halo of cold gas clouds. 
In principal, small cold gas clou ds from the IGM can be ac- 
creted onto the disks of galaxies (lBinneyiri977[ see also Sect. 1). 
However, the surrounding x-ray halo of massive galaxies ionises 
the c old gas before it can reach the disk (e.g. Nipoti & Binney, 
120071) . According to simulations of e.g. Dekel & Birnboim 
(2008), accretion of cold gas is not expected at galaxy masses 
above ~ 10 12 M . The dynamical mass of NGC 1167 is at least 
M dyn > 1.6 x 10 12 M Q (for r max = 65 kpc and v rot = 330 km s" 1 ). 
Hence, NGC 1167 may simply be too massive for the accretion 
of small gas clumps. 

Another possible explanation for the absence of an extended 
halo of cold gas clouds in NGC 1167 is the hypothesis that 
cold halos are caused, directly and indirectly, by star forma- 
tion. Because NGC 1167 does not form stars, galactic foun- 
tains from supernovae explosions do not blow any gas int o 
the halo. However, according to e.g. iMarinacci et al.1 (120101) . 
this galactic fountain gas would be needed to create Kelvin- 
Helmholtz instabilities in the hot corona such that (hot) gas can 
condense to clouds that are able to reach the disk. With the 
presence of supernovae, star-forming galaxies can accumulate 
enough cold gas to sustain th eir observed star formation rates 
(see also Hopki nset all l2008) while in early-type galaxies, like 
NGC 1167, the creation of instabilities in the halo is prevented 
(due to the absence of galactic fountains). Hence, no new cold 
gas complexes can reach the disk of NGC 1 167. 



6.3. Streaming motions, rotation curve and H 1 density 
distribution 

Streaming motions in disks are a well known phenomenon and 
observed in many galaxie s (see e.g. lKnapenill997l;|Pisano et all 
119981: iFathi et all 120061) . However, it is remarkable that the 
streaming motions in NGC 1167 are (to first order) symmet- 
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ric. The amplitude is relatively large 50 km s , depro- 
jected), at the high end o f what is detected i n a number of 
other early-type galaxies dNoordermeerl 120061) . Interestingly, 
most of these galaxies, like NGC 1167, also show clear signs 
of interaction and/or tidal distortions with some of them hav- 
ing bumps and wiggles i n the rotation curve as well (see also 
iNoordermeer et al.L l2QQ7h . 

The feature in the rotation curve of NGC 1167 occurs at 
the radius where the H i surface density drops significantly (see 
Fig. [2] and [5]). The drop in surface density results in the bump in 
the Hi scaling model rotation curve (Fig.fTUland Sect. 5.1), sug- 
gesting that a connection between large-scale streaming motions 
and Hi density distribution exists. 

H i scaling has previously been tried in the analysis of rota- 
tion curves and has even led to suggestions of a possi ble cou- 
pling o f dark m a tter wi t h the gas compo n ent (s e e. e.g.jBroeilsL 
1992[ ISwatersL 119991: iHoekstra et all 120011: iNoordermeerl 
20061) . Remarkably, in galaxies with bumps and wiggles, the 



rotation curves are fitted best by models with scaled H i disks 
(e.g. in UGC 2916, U GC 5253, UGC 6787 and UGC 1 1852; see 
iNoordermeeri 120061) . We therefore suspect that in many galax- 
ies a connection between H i density distributions and streaming 
motions in the disks may exist - because both are the result of 
an interaction. This may be the reason for the apparent success 
of H i scaling in describing such rotation curves. 



7. Summary 

We reported deep WSRT Hi observations of the early-type 
galaxy NGC 1 167. Our main conclusions are: 

(i) NGC 1 167 is a very massive early-type galaxy with a gi- 
ant Hi disk (£> ffl = 160 kpc). The Hi mass of 1.5 x 10 10 M 
contributes about 1% to the total mass, the Hi surface density 
(< 2 M Q pc~ 2 ) is too low for starformation to occur. 

(ii) The H i disk shows regular rotation for r < 65 kpc with a 
flat rotation curve of -340 km s" 1 amplitude. At r « 1.3 XR25 the 
rotation curve shows a bump of -50 km s" 1 on the approaching 
and receding side which is due to large-scale streaming motions 
in the disk. We argue that both the streaming motions and density 
features are due to interactions and that this explains the success 
of Hi scaling in fitting rotation curves. 

(iii) NGC 1 167 shows various signs of recent satellite accre- 
tion, interaction and merging activities as well as eight H 1 com- 
panions (with r < 350 kpc). No extended halo of cold gas clouds 
is detected with column densities above 6 x 10 18 cm -2 . This sug- 
gests that the H 1 disk is currently primarily formed through the 
accretion of companion galaxies which bring in large amounts 
of Hi to NGC 1167 and that the accretion of small cold gas 
clumps is not an important disk building mechanism in this mas- 
sive galaxy. 

In a future paper (Struve et al. in prep.) we will discuss 
whether our results, for instance the presence of streaming mo- 
tions in the disk, also have implications for the current phase 
of nuclear activity of the radio source B2 0258+35. The origin 
and onset of nuclear activity in galaxies is still not fully under- 
stood. Gas certainly plays an important role in the triggering of 
the black hole and NGC 1167 is a suitable object to study the 
interplay between the host galaxy and the radio source. 
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Appendix A: Channel maps Table B1: Rotation curve. 

Radius ["] Radius [kpc] velocity [km s _1 ] 
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Fig. A.l: Channel maps showing every third velocity channel of the data cube. Contour levels: -24, -12, -6, -3, -2 (turquoise), 2, 3, 
6, 12 and 24o- (black). 
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Fig. A.l: — Continued. 
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Fig. A.l: — Continued. 



C. Struve et al.: Cold gas in massive early-type galaxies: The case of NGC 1167, Online Material p 5 



16 
14 



16 
14 
12 
10 
35 08 



16' 
14 
12 
10 
35°08' 



■ 5116.39 . • 

© © 

to©. , w 

a v • 0: 

- 09 * " 9 i §> O 

- y), rt I i I a Al i I 


{f 



^ o - 

• $ - 


-"• 1 9 ' 1 1 vi',- 
5142.02 ^» y ' o J " 

- • ' V, * ° \ i 

■ °.'vio ■ 


5167.85? o • & % a 

: 4 » $ ' i^a£V ^ 
■ o m o; 

*. 'W\. ; / 
°° ft 

o * © ° % " 

'°- i ii^ +■ -i — i * J* ^ : 


" ■ lo ' — ^ — i — 1 ■ ' l j 

\ » @ 
6 ° ' ♦ °. « 

. o 

i . ■ A ' < 
• O-Cv fr> 

^ 9 

■ ' °9 o^l ■ 1 ■ ^1 ft> ^ 1 - 


$5218. 93rf 

• ° oo o » § ; ° . * 

o • * °' ' ' *° < 

A ° • 

• • ; • o o % 

o o n , 6 

V 

i i " i iO/ftii i • - 


5244.58° ° C'tO.V 

1 • # * f) 

^ • 

b 9 o * e 

o . • . ^ % - | 

: \ r V : 

" « 9 <? • ° « ° 



3 h 2 m s 50 s 40 s 30 s 20 s 50 s 40 s 30 s 

RA (J2000) RA (J2000) 



Fig. A.l: — Continued. 



